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I.  Introduction 

~  .,.V. 

Heat- transfer  by  convection  from  a  rotating  body  is  of  impc^^e 

A, 

in  the  therrr.al  analysis  of  rotating  components  of  various  types  oW«!^v:|^ 
chinerjr.  The  rotating  disk  is  one  of  a  number  of  geometrical  conf^^Js^^ 
tions  of  interest  because  many  practical  systems  can  be  idealized  jj^ 

terms  of  a  disk  rotating  in  an  infinite  environment  or  in  a  housing'^J. 

The  hydrodynamic  phenomena  associated  with  a  disk  rotating^^^jtjan' 
infinite  environment  have  been  investigated  theoretically  by  V. 

Cochran  (2),  Goldstein  (3),  and  Stuart  (li),  and  experimentally  by  li^Mper 
and  Regier  (5),  Smith  (6),  and  Gregory  and  Walker  (ii).  Heat  transyr  bj|||| 
convection  has  been  investigated  theoretically  by  Wagner  (7)  and 
and  Pohlhausen  (8)  for  laminar  flow  and  by  Kreith  and  Taylor  (9)  for 
turbulent  flow,  and  experimental  data  for  cooling  in  air  have  been  ob¬ 
tained  by  Young  (lO)  in  laminar  flow,  and  by  Cobb  and  Saunders  (ll)  in 
the  laminar  as  well  as  in  the  turbulent  flow  regime. 

For  a  disk  rotating  in  an  enclosure  Pantell  (12)  and  Jimbo  (13) 
have  measured  skin  friction  and  velocity  profiles  experimentally,  and 
Soo  (lU)  and  Schultz-Grunow  (l5)  have  investigated  the  boundary  layer 
and  flow  characteristics  theoretic^y;  neither  experimental  nor  theo¬ 
retical  work  on  heat  transfer  are  avjfcable  to  date. 


This  paper  contains  the  results  of  one  phase 
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program  designed  to  investigate  transfer  phenomena  in  rotating  systems. 

In  this  phase  of  the  work  mass  transfer  rates  from  a  disk  rotating  at  vari¬ 
ous  speeds  in  an  infinite  environment  in  laminar  and  turbulent  flow  were 
measured  and  compared  ^vith  the  results  of  an  analysis  which  is  also  appli¬ 
cable  to  heat  transfer.  The  effect  of  placing  a  stationaiy  surface  at 


various  distances  parallel  to  the  rotating  disk  surface  was  investigated 
experimentally  in  the  laminar  flow  regime.  Some  unusual  observations 
of  flow  patterns  resulting  from  transitional  "Goertler  Vortices"  and 
from  some  as  yet  unexplained  turbulent  vortex  phenomena  are  also  report¬ 
ed. 

II.  Experimental  Procedure 

During  the  initial  stages  of  developing  equipment  suitable  for 
studying  transfer  processes  from  a  rotating  disk,  an  electrically  neated 
test  section  was  considered.  The  design  of  this  test  section  involved 
the  use  of  a  guard  heater  at  the  nLm  of  the  rotating  disk  and  sliprings 
for  thermocouple  or  thermistor  leads  as  well  as  for  the  electrical  heat¬ 
ing  element.  The  initial  cost  of  constructing  this  type  of  test  section 
was  estimated  to  be  quite  high  and  the  reliability  of  measured  heat 
transfer  coefficients  appeared  to  be  uncertain.  In  order  to  circumvent 
the  experimental  difficulties  involved  in  taking  heat  transfer  and  sur¬ 
face  temperature  measurements,  it  was  decided  to  use  a  mass  transfer 
analogy,  similar  to  one  employed  by  Sogin  (l6).  This  decision  shifted 
the  experimental  difficulties  to  finding  a  suitable  material  which  could 
withstand  rotative  speeds  up  to  10,000  rpm,  but  eliminated  the  need  for 
taking  heat  flux  and  temperature  measurements  in  a  rotating  body,  and 
also  obviated  the  necessity  of  correcting  for  heat  leakages. 

The  arrangement  of  the  equipment  used  in  the  experimental  part 
of  this  study  is  shown  in  Fig.  1.  An  8  in.  diameter  aluminum  disk  was 
mounted  horizontally  on  a  Serval  "Superspeed"  centrifuge  which  co\fld 
be  adjusted  to  and  maintained  at  speeds  ranging  from  UOO  to  10,000  rpm. 


The  speed  was  measured  by  means  of  a  stroboscope.  The  temperature  of 
the  air  above  the  disk  was  determined  with  a  precision  thermometer.  The 
tests  were  carried  out  in  a  large  closed  room  to  eliminate  stray  air 
currents. 

Before  each  test  run  crystal  grade  naphthalene  molecular 

weight  123.2,  melting  point  79-81  C,  residue  after  ignition  less  than 
0.001/fc)  was  melted  in  a  covered  flask  in  an  electrically  heated  oven  at 
about  1201' .  The  molten  naphthalene  was  then  cast  into  the  undercut  and 
recessed  center  of  the  aluminum  disk  on  whose  rim  a  split  retaining  ring 
was  mounted.  The  disk  with  the  molten  naphthalene  was  then  cooled  in 
air  at  room  temperature.  After  the  naphthalene  had  solidified  the  re¬ 
tainer  ring  was  removed  and  the  surface  of  the  naphthalene  was  machined 
on  a  lathe  to  obtain  a  smooth  surface  raised  about  1/32  in.  above  the 
rim  of  the  disk. 

Prior  to  each  test  the  surface  of  the  naphthalene  was  carefully 
cleaned  of  loose  particles  and  the  aluminum  disk  with  the  cast  naphthalene 
coating  was  weighed  on  a  precision  balance  to  within  0.002  gm.  The  disk 
was  then  mounted  on  the  centrifuge  and  brought  up  to  speed.  The  time 
was  recorded  at  initial  weighing,  at  centrifuge  start,  at  attainment  of 
the  desired  speed,  at  centrifuge  power  off,  at  centrifuge  stop,  and  at 
reweighing. 

Since  the  measured  weight  difference  between  initial  and  final 
weighings  Included  losses  between  weighing  and  mounting,  between  dis¬ 
mounting  and  weighing,  and  during  the  start  up  and  slow  down  periods, 
appropriate  corrections  had  to  be  applied  to  obtain  the  difference  in 
weight  due  to  mass  transfer  at  the  test  speed.  The  losses  before  mount- 
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ing  and  after  dismounting  were  due  to  natural  convection  transfer.  Mass 
transfer  rates  during  natural  convection  were  therefore  measured  inde¬ 
pendently,  multiplied  by  the  total  time  elapsed  between  the  initial 
weighing  and  mounting  and  between  dismounting  and  the  final  weighing, 
and  the  resulting  quantity  subtracted  from  the  measured  weight  difference. 
To  correct  similarly  for  losses  during  speeding  up  and  slovfing  doim  of 
the  disk,  it  was  assumed  that  the  mass  transfer  rate  varied  linearly 
with  time  and  that  the  rates  of  mass  transfer  were  independent  of  the 
acceleration.  These  assumptions  were  checked  by  repeating  test  runs  at 
a  given  rotational  speed,  but  varying  their  total  duration.  Since  the 
results  were  reproducible  within  the  estimated  accuracy  of  the  essperi— 
mental  data,  the  method  of  correction  was  ccnsidered  satisfactory.  All 
test  runs  were  of  sufficient  duration  that  the  total  of  the  mass  loss 
corrections  necessary  for  any  of  them  was  less  than  1.5  pei*  cent  of  the 
measured  weight  difference. 

The  average  mass  transfer  coefficient  for  the  system  k^  is 
defined  by  the  equation 

where  m  is  the  total  rate  of  mass  transfer  in  Ib/hr 

is  the  ideal  gas  constant  of  the  vapor  in  ft  #/lb  R 
Tg  is  the  surface  temperature  of  the  solid  naphthalene  in  R 
P'yQ  is  the  vapor  pressure  of  the  naphthalene  at  Tg  in  #/sq. ft. 

A  is  the  surface  area  of  the  naphthalene  in  sq.  ft.  and 

kg  is  the  average  mass  transfer  coefficient  in  Ib/hr  sq.  ft.  (ib/ft^) 

*  The  concentration  of  naphthalene  vapor  in  the  air  above  the  plate  is 
so  small  that  no  correction  need  to  be  made  to  distinguish  the  process 
from  eqiiimolar  counterdiffusion. 
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In  Eq,  1  the  vapor  pressure  is  the  saturation  pressure  corresponding 
to  the  temperature  Tg  which  was  not  actually  measured.  However,  as 
shown  in  Ref.  17,  the  tenperature  measured  by  the  thermometer  differs 
from  the  disk  surface  temperature  only  by  two  small  corrective  terms. 

One  accounts  for  the  surface  temperature  depression  by  cooling  due  to 
sublimation.  Its  numerical  value  is  about  0.1  F.  The  other  term 
corrects  for  the  difference  between  the  recovery  temperature  indicated 
by  the  thermometer  and  that  of  the  cylinder  surface.  For  the  rotating 
disk  this  term  is  where is  the  recoveiy  factor. 

Assuming  that  the  recovery  factor  for  a  rotating  disk  is  of  the  same 
order  of  magnitude  as  for  a  flat  plate  at  the  same  velocity,  this  cor¬ 
rection  term  never  exceeds  1  degree  F.  The  omission  of  these  correction 
terms  does  therefore  not  introduce  an  appreciable  error  in  the  results. 
An  analysis  of  the  accuracy  of  the  experimental  restilts  in  terms  of  un¬ 
certainty  intervals  based  on  20  to  1  odds  (18),  indicates  that  the  value 
of  any  of  the  mass  transfer  coefficients  is  accurate  to  within  i  6  per 
cent  without  applying  a  surface  temperature  correction. 

The  physical  properties  used  in  the  reduction  of  the  experimental 
data,  including  the  diffusion  coefficient  of  naphthalene  vapor  in  air, 
the  specific  gravity  of  solid  naphthalene,  and  the  vapor  pressure  of 
solid  naphthalene,  were  taken  from  Table  2  of  Ref.  19  which  is  based  on 
the  best  available  data. 

In  the  tests  designed  to  study  the  effect  of  a  shroud  on  the 
transfer  process,  the  experimental  procedure  described  above  was  fol¬ 
lowed,  but  after  the  disk  was  mounted  a  glass  plate  was  placed  parallel 
to  and  at  a  predetermined  distance  above  the  disk  surface  on  an  adjust- 
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able  stand.  The  same  data  reduction  method  as  for  the  free  disk  was  used. 

Ill .  Evaluation  of  Average  Iviass  and  Heat  Transfer  Coefficients  for  a 
Disk  Rotating  in  an  Infinite  Environment. 

The  experimental  technique  used  in  this  investigation  as  well  as 
that  in  the  heat  transfer  measurements  of  Cobb  and  Saunders  (ll)  yielded 
average  coefficients  of  mass  and  heat  transfer.  To  compare  the  measured 
values  with  those  calculated  analytically  it  is  therefore  necessary  to 
evaluate  average  values  of  the  coefficients. 

As  shown  in  Appendix  I  as  well  as  in  Refs.  7  and  8,  the  theory 
predicts  that  in  the  laminar  flow  regime  mass  and  heat  transfer  co¬ 
efficients  are  uniform  over  the  entire  surface  of  a  rotating  disk. 

For  the  system  used  in  the  experiments,  i,  e.,  for  a  Schmidt  number 
of  2,14.3  the  approximate  analysis  presented  in  Appendix  I  yields  the 
equation 

Sh  =  ^7 (2) 

Equation  2  applies  also  to  heat  transfer  for  a  fluid  having  a  Prandtl 
number  of  2,h  if  the  Sherwood  number  is  replaced  by  the  Kusselt  number 
hcr/k. 

An  exact  solution  of  the  complete  Navier-Stokes  and  energy  equa¬ 
tions  for  the  rotating  disk  was  obtained  by  Millsaps  and  Pohlhausen  (8). 
For  cases  where  viscous  dissipation  can  be  neglected,  the  usual  situation 


in  laminar  flow,  their  results  can  be  reduced  to  the  equation 

A/u  -  -  C 


(3). 
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where  the  constant  C  depends  on  the  Schmidt  or  Prandtl  number  as  shown 
in  Table  I 

Table  I 

Pr  or  Sc  0.7ii  1.0  2.5  5.0  7.5  10 

C  in  Eq.  3  0.33  0.39  0.60  0.80  1.0  1.1 

As  shown  in  Ref.  1;,  laminar  flow  exists  at  Reynolds  numbers 
below  about  2  x  10^.  At  about  this  value  of  Reynolds  number  dis¬ 
turbances  in  the  laminar  boundary  layer  on  the  disk  are  amplified  and 
rapidly  cause  transition  to  turbulent  flow.  At  a  Reynold's  number 
larger  than  about  2.8  x  10^  the  flow  is  completely  turbulent.  On  a 
rotating  disk  having  a  Reynolds  number  larger  than  2.8  x  10^  there 
exists  therefore  a  laminar  inner  core  surrounded  by  a  ring  shaped 
transition  region  and  a  ring  shaped  turbiiLent  regime. 

In  the  turbulent  regime  mass-  and  heat  transfer  coefficients  in¬ 
crease  with  increasing  radius.  As  shown  in  Appendix  11^ a  modified  mass, 
heat,  and  momentum  transfer  analogy  predicts  local  values  of  Sherwood 
or  Nusselt  numbers  given  by  the  equations 


respectively,  where  the  symbols  are  defined  in  the  nomenclature. 

To  evaluate  average  values  of  the  mass  coefficients  the  local 
values  must  be  integrated  over  the  entire  disk  surface  according  to 


-li¬ 


the  equation 


where  r  is  the  radius  corresponding  to  the  critical  Eeynolds  number, 
c 

The  same  procedure  must  be  followed  to  evaluate  the  average  heat 
transfer  coefficient.  Since  the  expressions  for  and  h^^  derived  in 
Appendix  II  can  not  be  integrated  analytically,  the  averaging  was  done 
numerically  according  to 


where  ^r  was  taken  as  0-2  and 

N  =  (ro  “ 

In  the  averaging  procedure  the  transition  region  was  neglected 
and  it  was  assumed  that  turbulent  flow  sets  in  at  the  radial  distance 
corresponding  to  the  critical  Reynolds  n\miber  (see  Fig.  2).  The  average 
Sherwood  number  corresponding  to  the  average  value  of  the  mass  transfer 
coefficient  was  then  evaluated  according  to 

sjh 

A  similar  procedure  was  followed  in  the  evaluation  of  the  average 


Nusselt  numbers. 
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IV.  Presentation  and  Discussion  of  Results 

The  experimental  results  for  mass  transfer  from  the  disk  rotating 
in  an  infinite  environment  are  summarized  in  Table  II.  TJie  speed  was 
varied  from  about  500  to  10,000  rpm,  corresponding  to  a  rotational 
Reynolds  number  variation  from  about  3.5  x  10^  to  7  x  10^.  At  the 
highest  Reynolds  number  about  70  percent  of  the  total  disk  surface 
was  covered  by  a  turbulent  boimdary  layer. 

The  test  results  are  presented  graphically  in  Fig,  3  where  the 
average  Sherwood  number  is  plotted  as  a  function  of  the  rotational 
Reynolds  number.  Also  sho^vn  in  Fig.  3  are  the  results  of  theoretical 
analyses  for  the  laminar  flow  regime  and  the  results  of  the  analogy 
calculations  according  to  Eqs.  7  and  8.  In  addition  to  the  mass  trans¬ 
fer  results,  averaged  heat  transfer  data  obtained  by  Cobb  and  Saunders  (ll) 
with  a  disk  rotating  in  air  are  compared  with  the  results  of  the  analogy 
calculations  in  Fig.  3. 

In  the  laminar  regime  the  experimental  data  fall  slightly  below 
the  curve  predicted  by  the  analysis  in  Appendix  I,  but  slightly  above 
Millsaps  and  Pohlhausen's  theory  (8).  When  the  Reynolds  number  exceeds 
the  critical  value  and  part  of  the  disk  surface  is  covered  by  a  laminar 
boundary  and  part  by  a  turbulent  boundary  layer,  the  experimental  mass 
transfer  data  fall  slightly  below  the  results  predicted  by  the  analogy 
if  a  critical  Reynolds  number  of  2  x  10^  is  used,  but  above  the  calcu¬ 
lated  curve  if  2.5  x  10^  is  taken  to  be  the  critical  Reynolds  nmber. 

In  view  of  the  numerous  assumptions  made  in  the  analysis  the  agreement 
between  theory  and  experiment  is  quite  satisfactory  for  a  Schmidt  number 
of  2.U. 
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SUMMRY  OF  RESULTS  FOR 
MASS  TRANSFER  FROM  DISK  ROTATING 
IN  AN  INFINITE  ENVIRONMENT 


Test 

Run  No. 

Speed 

(RPM) 

Temperature 
of  Air  (®F) 

Rate  of  Mass 
Transfer 
(GM/Min.) 

Average 

Sherwood 

Number 

1 

2.900 

79.7 

0.02546 

293 

2 

4.500 

77 

0.04074 

545 

3 

5.500 

77 

0.05776 

773 

4 

2.000 

77.4 

0.01762 

231 

5 

7,000 

77.4 

0.08817 

1,154 

6 

7,200 

76.6 

0.08520 

1,216 

7 

1,000 

75.9 

0.01137 

162 

8 

2,000 

72.9 

0.01421 

240 

9 

3,000 

72.9 

0.01812 

305 

10 

6.520 

74.9 

0.06238 

937 

11 

500 

69.8 

0.00568 

114 

12 

1,200 

72.1 

0,0101 

179 

13 

3,600 

79.6 

0.03364 

388 

14 

3.300 

77.9 

0.02931 

373 

15 

8,300 

78.0 

0.10949 

1,397 

16* 

10,300 

77.2 

0.I5IU 

2,005 

17 

700 

72.8 

O.OO852 

144 

18 

1,500 

70.8 

0.01065 

202 

19 

1,730 

70.9 

0.01169 

220 

20 

2.250 

76.7 

0,01896 

258 

21 

2,600 

77.7 

0.02144 

275 

*In  this  test  a  small  piece  of  solid  naphthalene  chipped  off  near 
estimated  that  the  actual  mass  transfer  rate  is  about  6  per  cent  ! 
value • 


rim  and  it  was 
than  the  measured 


The  experimental  heat  transfer  data,  that  is  the  average  Nusselt 
numbers  for  the  disk,  agree  in  the  laminar  regime  within  5  percent  v/ith 
the  results  predicted  theoretically  in  Refs.  ?  and  8.  At  turbulent 
Reynolds  niimbers  the  heat  transfer  data  fall  as  much  as  10  per  cent 
below  the  result  predicted  on  the  basis  of  the  analogy  calculations 
for  a  Prandtl  number  of  O,?!;  and  -A  of  unity.  One  possible  reason 
for  this  deviation  is  that  the  analogy  calculations  take  no  account  of 
variations  in  physical  properties  which  are  insignificant  in  mass  trans¬ 
fer,  but  should  be  taken  into  account  in  heat  transfer.  Moreover,  since 
the  disk  with  which  the  heat  transfer  data  were  obtained  was  electrically 
heated,  the  surface  temperature  may  not  have  been  unifoimi.  It  is  of 
course  also  possible  that  the  assumptions  made  in  the  analogy  calcula¬ 
tions  could  be  responsible  for  the  deviation  between  theoiy  and  experi¬ 
ment. 

Some  observations  incidental  to  the  primary  objective  of  the 
research  are  worth  mentioning.  At  the  end  of  several  tests  in  the 
turbulent  Reynolds  number  range  the  imprints  of  certain  features  of 
the  flow  co-uld  be  observed  on  the  naphthalene  surface.  These  traces 
corroborate  the  observations  made  by  Gregory  and  Walker  (4)  means 
of  the  china-clay  evaporation  technique  in  the  transition  regime. 
According  to  these  authors  a  vortex  system  is  generated  in  the  region 
of  flow  upstream  of  transition,  and  between  two  critical  radii  vortices 
occur  which  are  stationary  relative  to  the  moving  surface.  Also  the 
impilnts  left  on  the  naphthalene  surface  indicate  the  presence  of  fairly 
large  scale  vortex  patterns  locked  to  the  surface  in  the  transition 
regime.  However,  whereas  in  the  British  experiments  (U)  no  traces  of 
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the  vortex  pattern  -were  visible  in  the  turbulent  regime,  the  imprints 
on  the  naphthalene  suggest  that  also  in  the  fully  developed  turbiilent 
regime  vortices  remain  attached  to  the  disk  surface  and  penetrate  the 
so  called  laminar  sublayer  in  a  more  or  less  regular  manner.  The  scale 
of  the  vortices  in  the  turbulent  regime  appears  to  be  much  smaller  than 
in  the  transition  regime.  To  verify  the  existence  of  stationaiy  vortices 
in  the  turbulent  zone  the  tests  should  be  repeated  with  a  disk  whose  sur¬ 
face  is  covered  with  china-clay,  but  instead  of  looking  for  traces  merely 
after  completion  of  a  test,  the  streak  pattern  should  be  recorded  photo¬ 
graphically  throughout  the  test.  It  is  planned  to  perform  such  tests 
in  the  future. 

The  effect  of  a  shroud  on  the  mass  transfer  process  from  a  rotating 
disk  was  studied  e:jq)erimentally  in  the  laminar  flow  regime  at  rotational 
speeds  of  1,000  and  2,000  rpm,  corresponding  to  rotational  Reynolds 
numbers  of  about  7  x  10^  and  l.!i  x  10^.  The  test  results  are  simimarized 
in  Table  III  and  presented  graphically  in  Rig.  U  where  the  ratio  of  the 


average  Sherwood  number  for  the  shrouded  disk  to  the  Sherwood  number  for 
the  free  disk  at  the  same  rotational  Reynolds  number  is  plotted 

versus  the  ratio  of  the  distance  between  the  disk  and  the  shroud  surface 


and  the  disk  diameter,  y/d. 

An  inspection  of  the  results  shows  that  within  the  range  of 
Reynolds  number  covered  in  the  experiments  the  effect  of  the  Reynolds 


nvmiber  on  the  Sherwood  number  ratio  is  less  than  the  scatter  of  the  data. 
The  presence  of  a  shrouding  surface  decreases  the  rate  of  mass  transfer. 

At  a  clearance  to  diameter  ratio  of  0.25  the  mass  transfer  for  the  shroud¬ 
ed  disk  is  only  about  10  percent  less  than  that  for  a  comparable  free  disk 
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TABLE  III 

SUMMAKf  OF  RESULTS  FOR  MASS  TRANSFER 
FROM  A  DISK  ROTATING  IN  A  SHROUD 


Test  Run 

No. 

Speed 

(RPM) 

Dimensionless 
Shrouding  Distance 
(y/d) 

Sherwood  No. 

Ratio 

(Shy/Sh^) 

a 

1,000 

0.11i3 

0.871 

Id 

1,000 

0.059 

0.823 

c 

1,000 

0.129 

0.878 

d 

1,000 

0.031 

0.717 

e 

1,000 

0.09li 

0.872 

f 

1,000 

0.250 

0.896 

h 

2,000 

0.079 

0.893 

i 

2,000 

0.031 

0.75U 

3 

2,000 

0.035 

0.776 

k 

2,000 

0.023 

0.7U 

1 

2,000 

0.016 

0.651 

m 

2,000 

0.012 

0.585 

n 

2,000 

0.250 

0.926 
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and  drops  to  80  percent  of  the  free  disk  mass  transfer  rate  at  a  clear¬ 
ance  ratio  of  0.05«  As  the  clearance  is  further  decreased  the  mass 
transfer  rate  decreases  rapidly. 

These  observations  are  in  qualitative  agreement  with  analyses  of 
the  boundaiy  characteristics  in  this  sort  of  system  (lii,  15)»  Wien  the 
fixed  boundary  is  sufficiently  far  away  from  the  surface  of  the  rotating 
disk  there  will  be  two  distinct  boundary  layers,  one  on  the  stationary, 
the  other  one  on  the  rotating  surface,  leaving  between  them  a  "core"  of 
finite  width  which  rotates  at  approximately  one  half  the  velocity  of  the 
disk.  In  this  regime  the  presence  of  the  stationary  surface  v/ould  not 
be  expected  to  affect  the  heat  or  mass  transfer  appreciably.  However, 
as  the  clearance  is  decreased  a  point  is  reached  where  the  two  boundary 
layers  meet.  Any  further  decrease  in  clearance  will  impede  the  circula¬ 
tion  appreciably  and  eventually  cause  choking  of  the  central  core.  In 
this  regime  the  rate  of  heat  or  mass  transfer  would  be  e3q)ected  to  drop 
off  rapidly  with  decreasing  clearance.  The  experimental  results  corrob¬ 
orate  these  phenomena  qualitatively,  but  additional  studies  of  the  flow 
and  boundaiy  layer  characteristics  are  necessary  before  the  various  flow 
regimes  can  be  delineated  quantitatively. 
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2,  In  the  laminar  regime  of  rotational  Reynolds  numbers,  a 
stationary  surface  placed  parallel  to  the  surface  of  a  rotating  disk 
reduces  the  rate  of  mass  transfer,  but  the  reduction  is  less  than 

20  percent  of  the  transfer  rate  for  a  free  disk  at  clearance  to 
diameter  ratios  larger  than  0.05- 

3.  Traces  left  by  the  air  f leaving  over  a  disk  rotating  in  an 
environment  suggest  that  stationary  vortices  occur  in  the 

turbulent  as  well  as  in  the  transition  regime. 
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NOMENCLATURE 


The  following  nomenclature  is  used  in  this  papers 

A  ■  heat  or  mass  transfer  area^  sq.  ft* 

C  -  local  drag  coefficient  at  r 
Dr 

-  average  drag  coefficient,  — 2^  rdr 

^  yr  r^2 

Cp  -  specific  heat  at  constant  pressure,  Btu/lb  F 
d  *  diameter  of  disk,  ft* 

“  diffusivity  of  vapor  in  mass  transfer,  sq.  ft ./hr. 

■  conversion  factor,  32*2  x  (3600)®  ft*  lb/#  hr.® 
hg  local  convection  heat  transfer  coefficient  at  r,  Btu/hr*  sq*"ft*  F 


M 

m 

Pv 

P 

Q 

r 


o 

Rv 


■  average  convection  heat  transfer  coefficient, 

J'^o  hg  2yf  rdr,  Btu/hr*  sq.“ft*  F 

»  thermal  conductivity,  Btu/hr*  ft*  F 

■  local  im  ss  transfer  coefficient  at  r,  Ib/hr*  sq*  “ft*  (lb/cu.”ft*) 

■  average  mass  transfer  coefficient 

*  ■  \ —  r^o  kc  2lf  rdr,  Ib/hr*  sq.-ft*  (lb/cu*"ft*) 

7/  *‘o‘^  -b 

-  mass-flux  paraiMter  (a/A)^  Ib/hr*  sq.ft. 

-  rate  of  mass  transfer,  Ib/hr* 

•  partial  pressure  of  vapor  component  A,  #/sq*ft. 

»  total  pressure,  #/sq.ft, 

■  rate  of  heat  transfer  by  convection,  Btu/hr. 

■  heat  flux  parameter  (q/A)8/cp/3  ,  ft*  F/hr. 

-  radial  distance,  ft* 

■  critical  radius  at  which  transition  from  laminar  to  turbulent 

flow  begins,  ft. 

■  outside  radius  of  disk,  ft* 

*  ideal  gas  constant  for  naphthalene  vapor,  ft*  #/lb.  R 


Nomenclature  (contlnixed) 
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S  - 
T  - 

V  (r) 

^  * 

u 

7 

oC 

S 


i  shear  velocity,  ft. /hr* 
temperatiire,  F 

«  net  velocity  at  r  relative  to  disk  velocity  at  r,  ft, /hr. 

radial  velocity  ccmponent,  ft ,/hr, 

peripheral  velocii^  component,  ft ./hr* 

distance  Arcmi  disk  surface,  ft, 

thermal  diffusivity,  k/^c^,  sq.“ft,/hr. 

eddy  diffusivity  for  momentum,  sq,~ft,/hr, 

eddy  diffusivity  for  heat,  sq.“ft,/hr, 

eddy  diffusivity  for  mass,  sq.^ft./hr, 

hydrodynamic  boundary  layer  thickness,  ft. 

thermal  boundary  layer  thickness,  ft,  ' 

mass  transfer  boundary  l^r  thickness,  ft, 

kinematic  viscosity,  ^  aq.-ft,/hr. 


6J 

r 

7* 

7 

T 

7* 

Nu 

IE 

Sc 

/^H 


rotational  speed,  rad./hr, 
fluid  density,  Ib/cu.  ft, 
fluid  shear  stress,  #/sq,“ft. 
fluid  viscosity,  #  hr./sq.-ft, 

nondimensional  velocity  in  the  boundary  layer,  v/S 
distance  from  disk  surface,  ft, 

clearance  between  disk  surface  and  stationary  shroud  surface,  ft. 
nondimensional  distance  from  the  surface,  yS/>> 
local  Nusselt  number,  h^  r/k 
average  Nusselt  number,  r^^ 

Schmidt  number,  ^ 

ratio  of  eddy  dlffusivities  of  mass  and  momentum 
ratio  of  eddy  dlffusivities  of  heat  and  momentum 


fw/Gi 
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Nomenclature  (continued) 

Sh  =  local  Sherwood  number,  k^r/D^ 

Sh  =  average  Sherwood  number,  k^rQ/l)^ 

Shy  =  average  Sherwood  nvunber  for  a  shrouded  disk  with  clearance  y 
St  =  local  Stanton  number,  h^/cO  r/c, 

St  =  average  Stanton  number  h^/cdr^J^c^ 

Re  =  local  rotational  Reynolds  number,  cDv^/-i) 

RSq  =  rotational  Reynolds  number,  6? rQ^i^ 

Pr  =  Prandtl  number,  T>  lot 


Subscripts 

H  -  heat  transfer 

1  -  laminar 

M  -  momentum  transfer 

m  -  mass  transfer 

o  -  outside  radius  of  disk 

s  -  surface  of  disk 

t  -  turbulent 

ao  -  far  away  from  disk 

b  -  buffer  layer 
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Appendix  I 

Approximate  Calculations  of  Rates  of  Heat  and  Mass  Transfer  from  a  Disk 
Rotating  in  an  Infinite  Environment  in  Laminar  Flow. 


In  cylindrical  coordinates  the  steady  state  equation  for  mass 
transfer  of  component  A  by  molecular  diffusion  through  an  inert  gas  in 
laminar  motion  is 

n  y.  J.  ^  ^ ^  _i!£.  ss  Q  1 

■where  is  the  diffusi-vity  of  component  A  in  ft^hr  and  c  is  the  con¬ 
centration  of  component  A  in  the  inert  gas  in  Ib/cu.  ft. 

The  corresponding  equation  for  heat  transfer  "to  or  from  a  ro'ta- 
ting  plate  is 

_  -\-r  “t-T- 

O  2 


r  dr 


if  -Tiscous  dissipation  is  neglected  and  -the  fluid  properties  are  unifom. 
Assuming  that  the  temperature  or  the  concentration  at  the  surface  of  the 
ro'bating  plate  is  uniform  and  that  the  -temperature  or  the  concen'bration 
in  the  fluid  far  away  from  the  plate  is  zero,  the  boundary  conditions 
pertaining  to  the  system  are 


at  y  =  0:  or  c  =  Cg 

at  y  =  oo  T  -  0  or  c  =  0 

Von  Karman  (l)  and  Millsaps  and  Pohlhausen  (8)  have  shown  that  the 
boundary  layer  thickness  S  is  constant  o-^er  the  entire  plate,  provided 
the  small  influence  of  flow  toward  the  plate  is  neglected;  and,  that  the 
radial  flow  velocity^!,  at  a  given  dis'bance  y  from  the  disk  surface  is 
proportional  to  r.  Consequently  the  total  mass  rate  of  fluid  flowing 
outward  in  the  radial  direction  at  any  distance  r  is  proportional  to  'the 
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area^  r^.  Thus  a  concentration  or  tenperature  distribution  dependent 
on  y,  but  independent  of  r  may  be  assmed  (7).  Equations  1  and  2  re¬ 
duce  then  to  the  ordinary  differential  equations 

and 

respectively.  Equations 3  and  U  yield  upon  integration, 

and 

Another  integration  gives 


«0  ^ 

T  exp  cfy^Cj  cfy 


The  integration  constant  C  can  be  evaluated  from  the  boundary  condition 
at  y  «  0  while  the  boundary  condition  at  y  »  oo  is  automatically  satis¬ 
fied  by  the  upper  limit  of  the  first  integral  in  Eqs.  7  and  8.  The 
concentration  or  the  temperature  at  the  plate  surface  are  according  to 
Eqs.  7  and  8  respectively 


■7;=fex/>C;^fe/^<^y^Qe/y 
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On  division  of  Eq.  5  by  Eq.  9  or  Eq.  6  by  Eq.  10  the  integration  con¬ 
stant  C  drops  out  and 

e/c  _  ^ 

or  ® 

y 

JT  _  -r 

For  naphthalene  vapor  diffusing  from  the  solid  naphthalene  - 
air  interface  at  the  disk  surface  into  the  surrounding  air  the  con¬ 
centration  in  the  gas  phase  is  related  to  the  partial  pressure  p^ 
by  the  equation 

<===’P^/^yT  13 


where  p^  is  the  partial  pressure  of  the  naphthalene  vapor  in  the  air 
in  psf,  R^is  the  ideal  gas  constant  in  ft  #/lb  E,  and  T  is  the  absolute 
teiqjerature  in  deg.  R.  The  mass  transfer  coefficient  and  the  heat 
transfer  coefficient  h^  are  given  by 


The  corresponding  dimensionless  numbers,  the  Sherwood  nvimber  and 
the  Nusselt  nxmiber  are 


I  -  ii 


A2/- 


Ar  __ 


.CO  y 


The  intep:rals  of  Eqs.  l6  and  17  can  be  evaluated  graphically 
with  the  aid  of  Von  Kannan's  solution  (l)  of  the  hydrodynamical  problem 
according  to  which 

y 

f or  y  4  ,  u  = 

^  c 

-  18 

j. 

and  for  y>^,  Uy  =  ^  0.^oS  Cl^cO)^  I9 

where  S*  ^  2,£S  20 

The  mean  temperature  during  the  experiments  was  F.  At 
atmospheric  pressure  and  74  F  the  kinematic  viscosity  of  air  ^  is 
0.60  ftVhr  (19)  and  the  vapor  diffusivity  is  0.29  ft^hr  (19),  so 
that  Sc  =  2,4.  With  this  value  for  Eq.  I6  becomes  upon  evaluation 
of  the  integral  in  the  denominator 

S/?  -  o.^y' 

Similarly,  for  Pr  =  2.4  the  Nusselt  number  is 

A/u  =  a  ^  22 

Wagner  (7)  evaluated  the  heat  transfer  coefficient  in  air 


(Pr  =  0.74)  and  found 


A/iz  =  0.3^  C -  0.34 


It  should  be  noted  that  both  heat  and  mass  transfer  coefficients  are 
miform  over  the  entire  surface  in  laminar  flov/. 
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Approximate  Calculation  of  Eates  of  Heat  and  Mass  Transfer  from  a  Disk 
Rotating  in  an  Infinite  Environment  in  Turbulent  Flow. 

The  analogy  between  heat-,  mass-,  and  momentxmi-  transfer  has  been 
used  successfully  to  predict  heat-  and  mass-  transfer  coefficients  from 
experimentally  measured  friction  coefficients.  Typical  analogy  solutions 
for  flow  inside  of  tubes  or  ducts  and  for  flow  over  plane  surfaces  are 
described  in  Refs.  21,  and  Kays  and  BJorklund  (25) 

obtained  recently  an  analogy  solution  for  a  circular  cylinder  rotating 
in  an  infinite  environment  and  found  that  the  solution  agreed  satis¬ 
factorily  with  experimental  data  for  air. 

'  The  analogy  method  is  based  on  the  assumption  that  the  transfer 
mechanisms  for  heat,  mass,  and  momentum  are  physically  similar,  or  can 
at  least  be  described  by  similar  equations  of  the  form 


for  the  momentum  transfer  and  resulting  shearing  stress,  by 


for  the  rate  of  heat  transfer,  and  by 


for  the  rate  of  mass  transfer,  where  ,  and  are  dependent 

on  the  amount  and  kind  of  turbulent  mixing  at  a  point.  To  complete  the 
analogy  one  assumes  that  and  where  and 

must  be  determined  theoretically  or  experimentally  inasmuch  as  the 
equation  of  motions  for  a  fluid  contains  terms  which  are  not  present  in 
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the  diffusion  equation.  In  the  vicinity  of  the  surface  the  shear  stress, 
the  heat  flux,  and  the  mass  flux  are  nearly  constant  and  if  the  fluid 
properties  are  uniform,  Eqs.  1,  2,  and  3  can  be  written  with  reference 
to  conditions  at  the  surface  as 

•if  1. 


The  application  of  the  analogjr  method  to  heat  and  mass  transfer 
from  a  rotating  disk  is  greatly  simplified  if  one  assumes  that  the  flow 
in  the  immediate  vicinity  of  the  disk  surface  is  similar  to  the  flow 
near  the  wall  of  a  pipe.  This  assumption  can  be  partly  Justified  on 
the  basis  of  available  data  on  the  drag  coefficients  in  the  turbulent 
regime  (5)  which  are  well  correlated  by  the  equation 

Equation  7  is  of  the  same  foim  as  the  well  known  Kaim^-JJ ikuradse 
resistance  equation  for  pipe  flow  and  the  coefficient  i;,07  corresponds 
to  a  magnitude  of  0,h  for  the  universal  non-dimensional  profile  constant 
for  turbulent  flow  near  a  surface.  Moreover,  measurements  by  Gregory 
and  Walker  (i;)  and  Theodorsen  and  Regier  (5)  indicate  that  the  velocity 
profile  near  a.  disk  is  logarithmic  outside  the  so  called  laminar  sub¬ 
layer  and  that  the  radial  velocity  component  is  of  a  smaller  order  of 
magnitude  than  the  tangential  velocity  component. 

Although  boundary  layers  on  rotating  disks  are  three  dimensional 
and  the  differential  equations  for  energy  and  momentum  as  well  as  the 
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boundary  conditions  are  dissimilar,  for  the  purpose  of  the  following 
approximate  calculations  a  conventional  analogy  approximation  was  made. 
The  flow  field  was  divided  into  a  socalled  laminar  sublayer  extending 
from  y"^::  0  to  y’*‘=  5>  a  bufferlayer  from  y'*’=  5  to  y'*‘=  30,  and  a 
turbulent  region  beyond  -  30.  The  upper  limit  of  the  buffer  layer 
thickness  was  estimated  by  comparing  the  temperature  distribution  pre¬ 
dicted  on  the  basis  of  strict  analogy  calculations  (25)  with  the  tem¬ 
perature  profile  measured  by  Cobb  and  Saunders  (ll).  The  solid  line 
in  Fig. S  represents  the  calculated  temperature  profile  while  the 
dashed  line  was  faired  through  esperimental  data  taken  from  Ref.  11. 

The  experimental  data  follow  the  predicted  curve  quite  well  up  to  a 
value  of  y  of  about  30,  but  beyond  this  point  the  erqperimental  tem¬ 
perature  gradient  becomes  less  than  the  calculated  value.  On  the  basis 
of  similar  observations  in  experiments  on  heat  transfer  from  a  rotating 
cylinder  Kays  and  Bjorklund  (25)  suggested  that  in  the  region  beyond 
y  =  20  the  transfer  mechanism  could  be  approximated  by  the  simple 
Reynolds  analogy.  Although  the  two  systems  are  dissimilar  in  mar^ 
respects,  the  Reynolds  analogy  was  found  to  be  a  satisfactory  approxi¬ 
mation  for  heat-  and  mass-  transfer  from  a  rotating  disk  beyond  y"*’=  30. 
Although  the  following  derivation  is  for  mass  transfer,  it  applies 
equally  well  to  heat  transfer  since  Eqs.  2  and  3  are  analogous. 

In  the  laminar  sublayer  y"*" ,  =  0  and  integration  of 

Eq,  6  gives  therefore 

AC^  =SMSc/S 
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In  the  buffer  layer  between  y  +  =  5  and  y+  =  30  the  von  Karman- 
Nikuradse  profile  (21)  is 

- 3.0S 9 

which  gives  a  dimensionless  velocity  gradient 

substituting  5/y+  for  the  velocity  gradient  in  Eq.  and  solving  for 
yields 
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Eq.  6  can  then  be  written  as 


so  that 
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Integration  between  y+=  5  and  30  yields  the  concentration  potential 
drop  in  the  buffer  layer  ,  or 

A/^  ^  1  10 


In  the  flow  regime  outside  of  y+=  30  it  will  be  assumed  that  the  inten¬ 
sity  of  transfer  processes  by  eddy  diffusion  is  much  greater  than  by 
molecular  diffusion  and  that  the  shear  stress  and  mass  flux  are  the  same 
function  of  y  so  that 

_ ^ 

^  a/  c: 

Integration  between  y"^  =  30  and  oo  yields  ,  the  concentration 
potential  at  the  outer  edge  of  the  buffer  layer,  or 
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Adding  the  concentration  potentials  and  noting  that  eOr/s  -  gives 


S/i  _  _  _ /^>7  y^£>/-/2 _  13 

For  the  calculations  in  this  paper,  the  experimental  results  of  Bhanak  (26) 
were  used  to  evaluate  /I^and  the  local  drag  coefficients  were  taken  from 


Ref.  5. 


Similarly,  the  analogy  gives  for  heat  transfer 


The  selection  of  an  appropriate  value  of/<7jj  for  the  numerical  evaluation 
of  Eq,  14  was  difficult  because  even  for  a  pipe  or  a  duct  available  experi¬ 
mental  results  are  inconclusive  (27,  28,  29,  30).  Experimental  measurements 
ofy^jj  taken  with  air  flowing  in  rectangular  channels  are  summarized  in  Ref, 

27,  They  indicate  that  the  space  average  value  increases  with  Rey- 

n 

nolds  number  from  about  1,1  at  a  Reynolds  number  of  2  x  10  to  1.4  at  a  Rey- 

<7 

nolds  number  of  9  x  10  but  varies  also  with  the  distance  from  the  surface 
at  any  given  Reynolds  number  (27,  28).  Data  obtained  with  liquid  metals  (29) 
give  values  of  range  from  0,7  to  0.9. 

In  the  past  a  value  of  unity  for^^  (20,22,25)  has  been  used  with  good 
success  for  heat  transfer  analogy  calculations.  Also  in  this  study  a  value 
of ®9^al  to  unity  was  found  to  yield  satisfactory  results,  giving  agree¬ 
ment  within  10  per  cent  between  the  Nusselt  numbers  calculated  from  Eq.  14  and 
the  experimental  results  of  fief.  11,  For  comparison  the  calculations  were  re¬ 
peated  with equal  to  1.4,  the  other  extreme.  This  caused  the  calculated 
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Nasselt  numbers  to  deviate  about  27  per  cent  from  the  measured  values* 

In  passing  it  is  interesting  to  note  that  semi-ea^jirical  methods 
deduced  originally  from  experimental  data  on  the  basis  of  a  thought  model 
implying  a  strictly  two-dimensional  structure  of  the  turbulent  boundary 
layer  have  in  this  work  and  in  Ref,  25  been  applied  to  systems  in  which 
the  boundary  layer  structure  must  be  three-dimensional  in  nature.  Since 
these  methods  are  apparently  applicable  to  cases  in  which  the  structure 
of  the  boundary  layer  can  by  no  stretch  of  the  imagination  be  considered 
two-dimensional,  it  seems  appropriate  to  question  the  hypothesis  of  a 
two-dimensional  turbulent  boundary  layer  in  the  other  cases.  The  pene¬ 
tration  of  the  so  called  laminar  sublayer  by  vortices  should  supply 
additional  food  for  thought. 
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